Polarizations in decays B u ^ — > VV and possible 
implications for R-parity violating supersymmetry 

Ya-Dong Yang 1 '* Rumin Wang 1,2 , G.R.Lu 1 ' 2 

1 Department of Physics, Henan Normal University, XinXiang, Henan 453007, P.R. China 
2 Institute of Particle Physics, Huazhong Normal University, Wuhan, Hubei 430070, P.R. China 

February 2, 2008 

Abstract 
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1 Introduction 



The study of exclusive nonleptonic weak decays of B mesons provides not only good opportu- 
nities to test the standard model (SM) but also powerful means to probe different new physics 
scenarios beyond the SM. Some important results from BABAR and Belle collaboration on the 
B decays to two light vector mesons B — > VV [TJI21IHllHElini0IElin] have caught considerable 
theoretical interests and have been extensively studied very recently [TUllTT|IT^lr^lT^lTo'l ll6|. 

In B — > VV decay modes, both longitudinal and transverse polarization states are possible. 
It has been known to us for many years that longitudinal polarization fraction is dominant in 
the SM when both the final vector mesons are light [TZ|. In the heavy quark limit, dynamics for 
the situation is analyzed clearly by Kagan ^H] and Grossman [TJ. Recent BABAR and Belle 
measurements of decay modes B — > p + p~, p°p + and p°K* + have confirmed this prediction as the 
final states are dominated by the longitudinal polarization. However, for the observed B — > <f)K* 
and p + K*° decays, the transverse polarization fractions are measured to be anomalous as large 
as 50%, which has posed a challenge to theoretical explanation. 

Recent studies have provided some possible resolutions to the anomaly. Kagan ^3] has 
shown it could be solved in the SM by increasing the nonfactorizable contributions of anni- 
hilation diagrams, which are formally power suppressed and depend on some poorly known 
nonperturbative QCD parameters. Later on, using a different factorization approach, Li and 
Mishima [H] have found that the annihilation contribution is not sufficient to lower the longi- 
tudinal fraction in B — > <pK* down to ~ 50%, although it could help to alleviate the anomaly. 
It is also shown that the anomaly might be due to large charming penguin contributions and 
final-state-interactions(FSI) by Colangelo et al. |T2] and Ladisa et al. However, Cheng 

et al. have found the FSI effects not able to fully account for the anomaly. Hou and Na- 
gashima [20] have given a model where the transverse descends from the transverse gluon from 
b — > sg*. However, it should be noted that <fi must couple to at least three gluons to neutralize 
color and conserve relevant quantum number. In this paper, we will explore the opportunity 
whether the RPV supersymmetry (SUSY) could provide a solution to the polarization anomaly 
and what kinds of constraints could be derived. 

The possible appearance of RPV couplings [21], which will result in lepton and baryon 
number conservation, has gained full attention in searching for SUSY [221 123- The rich 
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phenomenology implied by RPV couplings in B decays have been extensively discussed in 
|2H EE EE EH EH] • In this paper, we extend these studies to decays B — > VV. We use the 
QCD factorization(QCDF) framework [2E] for hadronic dynamics. We show that the polar- 
ization anomaly could be solved in the presence of RPV couplings, due to the appearance of 
q{V ± A)qb{V + A)s interactions. Compared with the matrix element of (K*\b(V — A)s\B), 
the axial parts in (K*\b(V + A)s\B) have opposite sign, hence induce polarization phenomena 
different from the SM. Moreover, using the recent experimental data for branching ratios and 
longitudinal-polarization fractions, we obtain limits on the relevant RPV couplings. 

This paper is organized as follows: In Sec. II, we calculate the B — > VV decay amplitudes 
which contain the SM part and RPV effects using the QCDF approach. In Sec. Ill, we list 
the theoretical input parameters used in our analysis. Section IV is devoted to the numerical 
results and discussions, we also display the allowed regions of the parameter space that satisfy 
all experimental data. Finally, our summary is presented in Sec. V. 



2 Amplitudes for B — > VV decays in QCD factorization 
approach 

Approaches for calculating amplitudes of B charmless nonleptonic decays always invoke fac- 
torization frameworks. Once factorization framework chosen, one can start from the AB = 1 
effective Hamiltonian of the underlying full electroweak (EW) theory at the renormalization 
scale /i ~ rrib, which can be obtained from the full EW theory by integrating out heavy particles 
with mass larger than using the renormalization group equation. The QCDF |2H] developed 
by Beneke, Buchalla, Neubert and Sachrajda is a powerful approach, which will be employed 
in this paper. Details of QCDF could be found in papers 

2.1 The decay amplitudes in the SM 

In the SM, the low energy effective Hamiltonian for the AB = 1 transition has the form [3U] 

nlfj = ^MCiO^+c^+^c^+c^^^ (i) 

V 2 i=3 
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with the effective operators given by 
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(2) 

Where v p = V* b V pq are the Cabibbo-Kobayashi-Maskawa (CKM) [HI] factors, Ci(fi) are the 
effective Wilson coefficients and could be found in Ref. [30] : a, (3 are the SU(3) color indices, 
{M2)v±a = qiY(± ± 7s)<?2, q = d,s and </ = it, d, s, c, 6. 

In order to calculate the decay amplitudes and branching ratios for B — > VV decays, we 
need the hadronic matrix element of the local four fermion operators 

V±A (hi) V-A \B), (3) 

where Ai, A2 are the helicities of the final-state vector mesons V\ and V 2 with four-momentum 
Pi and P2, respectively. In the rest frame of B system, since the B meson has spin zero, two 
vectors have the same helicity therefore three polarization states are possible, one longitudinal 
(L) and two transverse, corresponding to helicities A = and A = ± ( here Ai = A2 = A). We 
use X^ BVl ' V2 ^ to denote the factorizable amplitude for different chiral currents with the vector 
meson V 2 being factorized out. Under the naive factorization (NF) approach, we can express 
it as 

X (fl*,yy = { v 2 \(q 2llx (l -a 75 )g3)|0><MWl -7fc)fc)|£>, (4) 

here a = +1 or —1 corresponding to (q 2 q3)(v-A) or (q 2 q3)(v+A) current, respectively. 

We write the factorized matrix elements in term of the decay constant and form factors [22] 

(V(p,e*)\qi»q'\0) = f v m v e;, (5) 

(V(p,e*)\htf\B(p B )) = 2VBV } q2) e ^ £ * U PBP^ ( 6 ) 
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" m B + m v 
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q- 



-iq»(e*-p B )=^[A BV (q*)-A BV (q% (7) 



where ps(m^) is the four-momentum (mass) of the B meson, m^fej) and my 2 (e^) are the 
masses (polarization vectors) of the vector mesons V\ and V2, respectively, q = ps — p is the 
transferred momentum, and the form factors obey the following exact relations 
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One has the factorizable B — > V1V2 amplitude 
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mB + m Vl ~~ ''m B + m Vl 

We assume the V\ (V2) meson flying in the minus(plus) z-direction carrying the momentum p\ 
{P2), Using the sign convention e 0123 = —1, we get 
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V2) I - m B +m Vl v \" b V 2 > 

where ho for A = and h± for A = ±. 

The QCDF approach [29J allows us to compute the nonfactorizable corrections to the 
hadronic matrix elements (ViV2\Oi\B) in the heavy quark limit. The nonfactorizable corrections 
can be normalized to the factorizable amplitudes, so that they enter the effective parameters 
as a s corrections. They are calculated from the vertex corrections, hard spectator interactions, 
and QCD penguin-type contributions. 

In QCDF approach, light-cone distribution amplitudes (LCDAs) play an essential role. The 
LCDAs of the light vector meson are given in [22]. In general, the light-cone projection operator 
for vector mesons in momentum space can be divide into two parts 



M v = MY + Ml. 
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The longitudinal projector and the transverse projector are given by 
v f v m v (e*-n + ) v f^m v m v (e* ■ n + 
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Here we suppose the vector meson moving in the n_ direction, n± = (1,0, 0, ±1) is the light- 
cone vectors, u is the light-cone momentum fraction of the quark in the vector meson, fy and 
fy are vector and tensor decay constants, respectively, and E is the energy of the vector meson 
in the B rest system. In the main body of the paper we neglect power-suppressed higher-twist 
effects, i.e. we identify the meson momentum p' = En- and set e* ■ = 0. 

In the heavy quark limit, the light-cone projector for B meson can be expressed as (HH] 

ifBm B r 
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B 



with v = (1, 0, 0, 0) and the normalization conditions 
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where £ is the momentum fraction of the spectator quark in the B meson. 

The coefficients of the flavor operators a A which contain next-to-leading order coefficient 
and 0(a s ) hard scattering corrections can be written as follows: 
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where CV = (N^ — l)/(2Nc), Sj = mf/ml, and iVc = 3 is the number of colors. The superscript 
A denotes the polarization of the vector meson. 

In Eq.([15|). // (±1) contain the contributions from the vertex corrections, and given by 
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For hard spectator scattering contributions, we use the notation that V\ is the recoiled 
meson and V2 is the emitted meson, explicit calculations for fj T yield 
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with v — 1 — v. In Eq. (ll7|) . when the asymptotical form for the vector meson LCD As adopted, 
there will be infrared divergences in ffj. As in (HIEEE], we introduce a cutoff of order Aqct>/^6 
and take Aq C d = 0.5 GeV as our default value. 

The contributions of the QCD penguin-type diagrams can be described by the functions 
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with the function g(u, s) defined as 
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We have included the leading electro-weak penguin-type diagrams induced by the operators 0\ 
and 2 HE! 
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To calculate the coefficients a% in Eq. (|15jl . we have also taken into account the contributions 
of the dipole operator 8g , which are described by the functions 
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here we consider the higher-twist effects = uEnl + kj_ + ^^^+ in the projector of Eq. (|12|) . 
The G~ = in Eq. ((221 if we consider the Wandzura-Wilczek-type relations [SHj, but we get 
G+ ^ which is different from Ref. [1011211. 

With the coefficients in Eq. (|15|). we can obtain the decay amplitudes of the SM part A SM . 
B — > decay amplitudes are given in Appendix A. 

2.2 R-parity violating SUSY effects in the decays 

Flavor changing neutral current (FCNC) is forbidden at tree level in the SM. FCNC pro- 
cesses could happen at one loop level, but are suppressed by Glashow-Iliopoulos-Maiani (GIM) 
mechanism New physics effects could be comparable to the SM strength, so that penguin 
dominated rare B decays are sensitive to new physics. The RPV SUSY is an interesting scenario 
and its possible effects in B rare nonleptonic decays deserve our studies (221 12H] • 

The R-parity symmetry was first introduced by Farrar and Fayet which is assumed to 
forbid gauge invariant lepton and baryon number violating operators j2J. The R-parity of a 
particle field is given by R p = (— i) L + 2S '+ 3 - B ) where L and B are lepton and baryon numbers, 
and S is the spin. However, there is no deep theoretical motivation for imposing R-parity. The 
presence of RPV could give very rich phenomenology. Of course, it will get constraints from its 
phenomenology. The status of RPV SUSY and constraints on its parameters could be found 
in the recent reviews fZ*2\ 123] . 
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In the most general superpotential of the minimal supersymmetric Standard Model (MSSM), 
the RPV superpotential is given by jH] 

W$ = ^^,4^ + ^^ + ^^]^^, (23) 

where L and Q are the SU(2)-doublet lepton and quark superfields and E c , U c and D c are the 
singlet superfields, while i, j, and k are generation indices and c denotes a charge conjugate 
field. 

The A and A' couplings in Eq. (|2H|) break the lepton number, while the A" couplings break 
the baryon number conservation. There are 27 A'-type couplings and nine each of the A and 
A" couplings as \ij\k is antisymmetric in i and j, and is antisymmetric in j and k. The 
antisymmetry of the B- violating couplings, Af^w in the last two indices implies that there are 
no operators generating the b — > sss and b — > ddd transition. 

From Eq. (j2HJ), we can obtain the following four fermion effective Hamiltonian due to the 
exchanging of the sleptons: 

The four fermion Hamiltonian due to the exchange of the squarks can be written as 
A" A"* 
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Where P^, = ^-^,P R = = and /3 = H — The subscript for the currents 

(i/Ji,8 represent the current in the color singlet and octet, respectively. The coefficients t^ -4 /^ 
and ?7 -8 /^ are due to the running from the sfermion mass scale Vfip (100 GeV assumed) down to 
the mj scale. Since it is always assumed in phenomenology for numerical display that only one 
sfermion contributes one time, we neglect the mixing between the operators when we use the 
renormalization group equation (RGE) to run down to low scale. The ?fl for the relevant 
decay modes are written down in Appendix B. 



Compared with the operators in the Wfiy, there are new operators (q2Q3)v±A(bQi)v+A hi 
the T& 1 . We will use the (') denote the {q2(Jz)v±A(p<ii)v+A current contribution. In the NF 
approach, the factorizable amplitude can be expressed as 

X'WuV*) = (K 2 |(g 27At (l-a 75 )g 3 )|0)(V 1 |(67 At (l+75)gi)|5). (26) 

Taking the V\ (V2) meson flying in the minus (plus) z-direction and using the sign convention 
e 0123 = —1, we have 



if 



v 2 



X f(BV u V 2 ) 



2m 



[m% - m\ - m\){m B + m Vl )A? Vl « 2 ) ~ ^-^ Vl «, 



-if V2 m V2 



[m B + m Vl )A? Vl W 2 ) ± ^^V^(ml 2 , 



= /4. 



(27) 



Applying the QCDF approach and supposing V\ (V2) is the recoiled (emitted) meson, we 
obtain the vertex corrections f'j X (a) and hard spectator scattering corrections ///(a) as follows: 
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Since we are considering the leading effects of RPV, we need only evaluate the vertex 
corrections and the hard-spectator scattering. The R-parity violating contribution to the decay 
amplitudes A& can be found in Appendix C. 



2.3 The polarized fraction and branching ratio 

With the QCDF approach, we can get the total decay amplitude 

A x = Af M + Af. (29) 
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The expressions for the SM amplitude Af M and the RPV amplitude A\ are presented in 
Appendices A and C. From the amplitude in Eq. (|29|) . the branching ratio reads 

t~b\Pc 



Br(B -> VV) 



inm, 



i\A \ 2 + \A, 



(30) 



where tb is the lifetime of B meson, p c is the center of mass momentum, and given by 
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In order to compare the size of helicity amplitudes, we express the longitudinal polarization 
fraction 
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The ratios of Tl/T measure the relative strength of longitudinally polarizations in the decay. 



3 Input Parameters 

A. Wilson coefficients 

To proceed we use the next-to-leading Wilson coefficients calculated in the naive dimensional 
regularization (NDR) scheme and at mj scale |3*U] : 

Ci = 1.082, C 2 = -0.185, C 3 = 0.014, C 4 = -0.035, C 5 = 0.009, C 6 = -0.041, 
C 7 /a e = -0.002, C 8 /a e = 0.054, C 9 /a e = -1.292, C w /a e = 0.263, C 8g = -0.143. 

B. The CKM matrix element 

As for the CKM matrix elements, we will use the Wolfenstein parametrization [4*2] : 
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C. Masses and lifetime 

For quark masses, which appear in the penguin loop corrections with regard to the functions 
G x (s), we take 



m u = m d = m s = 0, m c = 1.47 GeV, m b = 4.8 GeV. 

To compute the branching ratio, the lifetime of B meson and masses of meson are also taken 
from [33] 

r Bu = 1.671 ps, m Bu = 5279 MeV, r Bi = 1.536 ps, m Bd = 5279 MeV , 
m K *± = 892 MeV, = 1019 MeV, m K ,o = 896 MeV, m p = 776 MeV. 

D. The LCDAs of the vector meson 

For the LCDAs of the vector meson, we use the asymptotic form ["""""] 



As for the two B meson wave functions in Eq. Ijlfljl . we consider only the $f (£) contribu- 
tion to the nonfactorizable corrections as done in the literature [SEllli]' since $f (£) is power 
suppressed. We adopt the moments of the $f (£) defined in Ref. [23 El] for our numerical 
evaluation: 



with \ B = 0.46 GeV ^j. The quantity A# parameterizes our ignorance about the B meson 
distribution amplitudes and thus brings considerable theoretical uncertainty. 
E. The decay constant and form factors 

The decay constant and form factors are nonperturbative parameters. They are available 
from the experimental data or estimated with theories, such as lattice calculations, QCD sum 
rules, etc. For the decay constant, we take the latest light-cone QCD sum rule results (LCSR) 
[3*2*] in our calculations: 



$[(*) = *v {x) = g wv = 
gW(x) = jj[l + (2*-l) 2 ]. 



6x(l — x), 



(34) 




(35) 



f K . = 217 MeV, f p = 205 MeV, fa = 231 MeV 
/£. = 156 MeV, = 147 MeV, f£ = 183 MeV, 
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and /s u(d) = 161 MeV. For the form factors involving the B —> K* and B — > p transition, we 
adopt the results given by jHZ| 

Af u(d)K \0) = 0.292, A? u(d)/< *(0) = 0.259, V B ^ K " (0) = 0.411, 
Af u(d)P (0) = 0.242, Af u(d)P (0) = 0.221, V B ^ P {0) = 0.323. 

4 Numerical results and Analysis 

In this section we will give our estimations in the SM and compare with the relevant experimen- 
tal data, then we show the RPV effects to the branching ratios and longitudinal polarization 
fractions. With the RPV effects we can do two things. First, it is possible to have small 
longitudinal polarization fractions for the pure penguin processes completely different from 
the SM predictions. Second, one can obtain more stringent bounds on the product combina- 
tions of RPV couplings if the measured values of the decay modes are consistent with the SM 
predictions. 

In the numerical analysis, we assume that only one sfermion contributes one time with a 
mass of 100 GeV. So for other values of the sfermion masses, the bounds on the couplings in 
this paper can be easily obtained by scaling them by factor f 2 = ( loocev ^ 2 - ^ e cons ider all 
possible sfermion contributions when we obtain the bounds on the couplings. 

4.1 B° -> (f>K*° and B+ -> cf)K* + 

B° — > <j)K*° and B + — > <pK* + decays are induced by the underlying quark transition b — > sss 
at one loop level. In the SM, these decays are called pure penguin processes. However, they 
could happen at tree level in SUSY models with RPV. Our estimations in the SM and recent 
data from the BABAR and Belle collaborations are presented in the Table I and II. 

From Tabled and II, we can see that the predictions from the two factorization framework 
are very similar. It is also noted that the latest LCSR results for factors are smaller than 
their pervious values j3H]- With the form factors, we find that the SM predictions for the 
branching ratios are smaller that the BABAR Ej and Belle [SI IE] measurements, and the 
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Table I: Branching ratios of B — > <pK* decay modes. 



Mode ^^^(xlO 6 ) Belle(xl0 6 ) Average(xl0 6 ) QCDF in SM(xl0 6 ) 

B+^(f)K* + 12.7l|g + l.ip 6.7l?]i?;S 9.5 ±1.7 4.60 

B° d ^<pK*° 9.2 ± 0.9 ± 0.5 [7] 10.0^1^0 9.4 ± 0.9 4.21 

Table II: The longitudinal polarization fractions of B — > <j)K* decay modes. 

Mode BABAR Belle Average QCDF in SM 

B+-^(pK* + 0.46 ± 0.12 ± 0.03 [I] 0.51 ± 0.08 ± 0.03 [H] 0.49 ± 0.07 0.861 

B* d <f)K*° 0.52 ± 0.05 ± 0.02 0.45 ± 0.05 ± 0.02 [B] 0.49 ± 0.04 0.861 



longitudinal polarizations are as large as ~ 90% in contrast to ~ 50% measured by BABAR 
and Belle. 

Now we turn to RPV effects which may give a possible solution to the polarization anomaly. 
Using the formula derived in previous sections, we can calculate RPV effects. The RPV effects 
to the branching ratios and longitudinal polarization fractions of the B — > <j)K* decays are 
displayed by curves in the Fig.l. We find that the A^A^ t erm could not enhance transverse 
polarization because the corresponding current has a structure S7 M (1 + 75)567^(1 — 75)5. Its 
matrix element is the same as the SM one, since (^sj^slO) = 0. However, the A^ 23 A^2 
term could provide a solution to the polarization anomaly because its current has the structure 
57^(1 — 75)367^(1 + 75)5. The right handed (bs) current will flip the signs of the axial parts 
of the matrix (K*\tr)^{\ — ^y 5 )s\B) appearing in the SM contributions. We also find that 
I A - 2 3 A 1 < 3.0 x 10 _3 z>|j by the branching ratios; however, there is a much stronger bound 
from the polarization ratios. 

Combining both the branching ratios and the polarization ratios which we get by QCDF 
(the same as follow), the resolution is obtained in a very narrow parameter interval |A^ 23 A^2 2 | £ 
[1.5 x 10~ 3 z>2j, 2.1 x 10~ 3 z>!j]. Fortunately, the parameter space is not ruled out yet, the existing 
upper limit is |A- 23 A-2 2 | < 2.3x 10~ 3 [221I211- We note that similar strength of the RPV couplings 
is also needed in the recent studies to solve the CP problem in B — > <pK [37| and the rj/ puzzle 
in B -> r]fK [22]. 



14 




Figure 1: The branching ratios and the longitudinal polarization fractions for B — > (pK* as functions 
of the RPV couplings A' A'*. The solid curves represent our theoretical results by QCDF. The dash 
lines are the predictions by NF for comparisons. The horizontal solid lines are the experimental data 
and the SM predictions with QCDF as labelled respectively. The horizontal dot lines represent the 
la error-bar of the measurements. (The same in Fig. 2 and 3). 

4.2 B+ -> pK* decays 

— > pK* decays are due to b — > uus or b — > dds transitions at the quark level. The 
B+ — > p + K*° decay is a pure penguin process in the SM. Its longitudinal polarization fraction 
were measured to be unexpected low ~ 0.5 very recently by Belle [H], which is inconsistent with 
the SM prediction. While the — * p°K* + has both tree and penguin amplitude, experimental 
measurements by BABAR £Q have shown the decay dominated by longitudinal polarization, 
which is consistent with the SM predictions. 

For comparison, our estimations in the SM and recent data from experiments are presented 
in the Table III and IV. 
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Table III: Branching ratios of B — > pK* decay modes. 



Mode 


BABAR (xlO 6 ) 


Belle(xl0 6 ) 


QCDF in SM(xl0 6 ) 


B + _> p °K*+ 
B + _> p+K *o 


10.6^ ± 2.4 P 


6.6 ±2.2 ±0.8 


2.58 
3.79 



Table IV: The longitudinal polarization fractions of B — > pK* decay modes. 



Mode 


BABAR 


Belle 


QCDF in SM 


B + -> p°K*+ 
B+ - p + K*° 


0.96^ ± 0.04 P 


0.50 ±0.19^ P 


0.906 
0.905 



From Table III, we can see the SM predictions for the branching ratios are smaller than the 
measurement by BABAR and Belle. As shown in Table IV, for the longitudinal polarization 
fraction in the tree dominant decay B + — > p°K* + , the SM prediction agrees with the BABAR 
measurement very well. However, for the pure penguin process B + — > p + K*°, the SM predicts 
dominant longitudinal polarization, which is in contrast to Belle measurement Pj. 

In the B+ -> p°K* + decay, since the quark content of p° is (uu — dd) / a/2, the decay could 
be induced by superpartners of both up-type and down-type fermions. For example, b —>■ dds 
could be induced by sneutrino, while b — > uus could be induced by slepton with the same 

A' A'* A' A'* — — — — 

A^ 13 A^ 2 product. We take !l V 12 and a V 12 contribute to b — > uus and b — > dds at the same 
time, so the effects of A^ 13 A^ 2 will be eliminated if taking = mo i . 

Our results for the RPV contributions to B+ — > pK* are summarized in Fig. 2. We can get a 
lot of information from this figure. We find that the longitudinal polarization in B+ — > p°K* + is 
sensitive to lA'^A'/^l and lA^A-^l arising from exchanging d and Ui between the quark currents 
n7 M (l + 7 5 )u ® &7 M (1 ±75)5 and d^^il ±75)0?® 67^(1 ±7s)s, respectively, and insensitive to the 
effect from exchanging z/j between the quark currents g?7 m (1 ± 75 )<i® 67^(1 —75)5. However, the 
pure penguin process B+ — > p + K*° is sensitive to the RPV couplings lA^A-'^l and (A^gA^-J 
which associate with operators dj^l ± 75) s <g> &7 M (1 ±75)0? and 6^(1 — 75)5 ® &7^(1 ±75)0?, 
respectively. 
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Figure 2: The branching ratios and longitudinal polarization fractions for B — > pK* as functions 
of the RPV couplings A' A'* and A" A"*. 

From Fig. 2, the polarization problem in B+ — » p°K* + could be solved by RPV effects. 
However, combining the constraints from polarization fraction and branching ratios of -B+ — > 
pK* , we obtain very narrow parameters spaces for relevant coupling constants. It implies that 
these couplings could be pinned down or ruled out by refined measurements in the very near 
future at BABAR and Belle. Our constraints are listed in Table V. 

For comparison, we have also listed existing bounds on these quadric coupling constant 
products. We can see that the product lA^A'/^l is severely constrained by double nucleon 
decay [211] • Our bounds on |A- n A^3 2 | is weaker than that by ArriK @B]- These couplings also 
contribute to B — > Kn decays. In an interesting study [27|, constraints on these couplings have 
been derived from the experiment measurement of branching ratios of B — > Kir. As shown in 
Table V, our bounds are lower than theirs |2Z|. However, in view of new data from BABAR 
and Belle, it would be interesting to investigate whether the RPV couplings could solve the 
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Table V: Bounds on the quadric coupling constant products. For comparison, the exiting 
bounds are listed in last column. 



Couplings 


Bounds 








Process 


Previous bounds 


^21-^31 


<5.1xl(T 3 i> 2 . 

— L% 

— 1 ,1. 


for 


A i21 A i31 <0 
A' A'* >0 


K 


p°A*+ 


8.2 x 10" 4 j22j 


I \" y* i 

l A i31 A i2ll 


<4.1xl0~ 3 u 2 

[4 Ox 1 _3 7/ 2 8 SxlO - 3 ?5 2 lanH<1 1 x10~ 3 ?7 2 


for 
for 


A i31 A i21 <0 

A'' A"* >0 
i31 i21^ u 




-> p°/l *+ 


1.0 x 10~ 2 |2l| 


1 \" \"* 1 
l A 131 A 12ll 


<1.8xlO" 2 J 2 

[1.6xl0- 2 d 2 ,3.8xl0- 2 d 2 ]and<4.7xl0" 3 (i 2 


for 
for 


A 131 A 121 >0 
A 131 A 121 <0 


Bt 


-> p°/C+ 


2 x lir 8 [23 


l A i23 A illl 


<1.5xl(r 3 i> 2 . 

[1.6xl0- 3 !> 2 l .,4.0xl0- 3 i> 2 l ]and<1.2xl0- 3 £; 2 l . 


for 
for 


A i23 A ill >0 
A *23 A ill <0 


Bt 


-> p+A*° 


2.1 x 1(T 3 j23 


l A ill A i32l 


<1.4xl0~ 3 i> 2 . 

— L% 

<4.1xl0" 3 i> 2 


for 
for 


A ill A i32 <0 
A ill A i32 >0 


K 


-> p+A*° 


4.7 x 10" 4 mi 


l A i31 A i2ll 


<2.5xlCT 3 « 2 

<6.3xl0- 3 u 2 

— i- 


for 
for 


A «l A £l <0 

A «i A £i>° 


Bt 


-> p + A*° 


1.0 x io~ 2 j23 



known irK puzzle |4*iH I5U] . 



4.3 B ^ pp decays 

Recent measurements of tree dominated vector- vector charmless modes —> p°p + and B% — > 
p + p~ by BABAR P3E] and Belle [2] have shown that the decays are dominated by longitudinal 
polarizations, which are just as the SM expectations. 

Our estimations in the SM and recent data from the BABAR and Belle collaborations are 
presented in Table VI and VII. We find that the SM predictions for both branching ratios and 
longitudinal polarization fractions agree with the BABAR and Belle measurements very well. 
Therefore, these decays could give strong constraints on relevant RPV couplings. 

The RPV contributions are presented in Fig. 3. From the figure, we can know that longitu- 
dinal polarization fractions for B — > pp are not sensitive to the RPV couplings since they are 
tree dominated decay in the SM. By the decay — > p°p + , we get the bounds on |A' 11 A^ 1 |: 

|\' \'* I c [6.3 xlO-^ 2 ,, S.OxlO- 3 ^] forA^A^X) , f>+ ^ n n + (oa\ 

l A ill A i3ll *= [8.6xl0- 4 £> 2 , 2.5xl0- 3 i> 2 ] for A' A'.* <0 °J D u ^ P P ■ W D ) 

L Lit 3 Li' ill lol 

The existing bounds on the products by B — > 7T7t decays are |A^ 11 A-3 1 | < 1.6 x 10 -2 [22j. The 
parameter spaces constrained by B£ — > p°p + are very narrow, which could be closed easily be 
future refined measurements if the tight bounds from branching ratio and polarization ratio 
do not overlap. This could not be done by the B — > PP, PV decays. We see again the rich 
phenomena of B — > VV decays and its power for bounding new physics. 
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Table VI: Branching ratios of B — > pp decay modes. 



Mode 


BABAR(xl0 6 ) 


Belle(xl0 6 ) 


Average(xl0 6 ) 


QCDF in SM(xl0 6 ) 


Bt - pV 
B° d - P + P - 


22.5±|;1 ± 5.8 [I] 
30 ± 4 ± 5 jU 


31.7 ± 7.1±l;f 


26.2 ±6.2 


15.30 
25.75 



Table VII: The longitudinal polarization fractions of B — > pp decay modes. 



Mode BABAR Belle Average QCDF in SM 

Bt -> pV 0.97lg:ol ± 0.04 [I] 0.948 ± 0.106 ± 0.021 ;2j 0.964 ± 0.056 0.941 
B° d -> p+p" 0.99 ± 0.03±g;g| jl] • • • 0.935 
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Figure 3: The branching ratios and longitudinal polarizations for B — > pp as a function of the 
RPV couplings A'A'* and A" A"*. 

We note that RPV coupling X'( 32 ^ii2 * s eliminated because the factorizable amplitudes 
X'( B ip°'p + ) = x'( BwP+ ' p °^ and cancelled each other. The effect of A^ 13 A^* 1 also is eliminated 
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for the same reason in B+ — > p°K* + . 

The decay 5° — > p + p~ give constraint \" 32 Ki2 e [9-3 x 10~ 4 s 2 , 1.1 x 10~ 3 s 2 ]. However, this 
parameter space is already ruled out in previous studies [221 [2Hj • Again, we see the standard 
model works well for tree dominant B decays. 

Through the numerical analysis in this section, we can conclude that the experimental value 
of the polarization anomaly in pure penguin decays can be solved with RPV effects; however, 
the solution parameter spaces are always very narrow. 

5 Conclusions 

Motivated by the polarization anomaly observed recently by BABAR and Belle, we have stud- 
ied B — > VV modes with QCD factorization approach both in the SM and in RPV SUSY 
theories. We have found a set of RPV couplings can give possible solution to the polarization 
anomaly. However, the windows of the RPV couplings intervals are found to be always very 
narrow. It implies that these couplings might be pinned down from the rich polarization phe- 
nomena in these decays. However, it also implies the window could be closed easily with refined 
measurements from BABAR and Belle in the near future. 

Since the hadronic dynamics for B nonleptonic decays are generally tangled with nonpertur- 
bative dynamics, we need to known how to separate perturbative and nonperturbative dynam- 
ics; of course, we also need to know the value for the parameters of nonperturbative dynamics 
to give reliable predictions based on electroweak theories. At the present stage, QCDF is a 
working scheme. Generally, we can believe that QCDF calculations for polarization fractions 
could be much more accurate than that for branching ratios, since many uncertainties could be 
cancelled in the fractions. Therefore the constraints from polarization measurements would be 
more well-founded than those from branching ratio measurements. 

Comparing our prediction with the recent experimental data, we have obtained bounds on 
the relevant products of RPV couplings. We find that many bounds are stronger than the 
existing limits fZ*2\ El] , which may be useful for further studying the RPV SUSY. 

In conclusion, we have shown that RPV SUSY could give possible solution to the polarization 
anomaly in pure penguin decays B — > <j)K* and B + — > p + K*° observed by BABAR and Belle. 
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Appendix 



A . The amplitudes in the SM 

A s x M (B^cPK*+) = -^V*V ts [a x + a x + a x 5 



A S X M (BZ - 

a s x m (b: - 



Af M (B° d 



<j)K*°) 
p°K*+) 



p + K 



*o i 



p + p~ 



V2 

~v; b v ts [a x 3 + 4 



a 9 + a 10 



A"* 



- a 7 + a Q + a 



10 



^ {KV us a> - V*V ts 3 -(a$ + al)]X^ + ^ 
+ K b V us a\ - V t tV ts (a x 4 + <)]X^°,^+)} 



Gj 



G 



[v: b V ud [a\X^^ 



-v; b v td 



^[V: b V ud a\ - V t tV td (a x 4 + a^X^'^. 



(37) 
(38) 

(39) 
(40) 

(41) 
(42) 



where we neglected the annihilation matrix element contributions. 



B . The Hamiltqnian for RPV 



2mi Li ^ — 2m^ 



p0ir+) = ^^^°{^^ 

-[(«7 M ^s)i(&7^)i " {ui^PRsUh^PRu)*]} 
-J^tr" Po [{d^PRd)i{h,P R s)i ~ (d^P R d) 8 (h,PRs)8] 

21 



v v* 

A„oi A„ q 



2V2m 



idYP R d)s(b^P L s) 8 V~ 



^f//(^ + ^P + ^°) = -f^V-^Kd^PRsU h,P R d)i - {d^P R s) 8 {h,P R d) 8 } 

111 

+ ^^V- 8/ "°(dl"PLs)s(h,P R d)s + ^^V- 8/ "°(d^P R sUh,P L d) 8 , 

\ll \//* 

^?//(^ + -pV) = ^yfP" 4/A) {[{u^PRuUh.Pitd), - (urP R u)s(h,P R d)s] 

-[(uYP R dUb % ,P R u)i - (u^PndUh^Pnu)* 

+ ( (u^PrvXJb^.P^ - -^k^ 



(uYP L u)s(b^P R d)s - % f (drP L d) 8 (b 7ll P R d) 8 



^^{drp R d) 8 {p %l p L d) 8 ^ v -w°, 



2m 



C . The amplitudes for RPV 

J$(B+ ^ <f>K* + ) = 
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